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We consider a phantom chain model of polymer with internal friction in a har¬ 
monic confinement and extend it to take care of effects of solvent quality following a 
mean field approach where an exponent n is introduced. The model termed as “Sol¬ 
vent Dependent Compacted Rouse with Internal Friction (SDCRIF)” is then used 
to calculate the reconfiguration time of a chain that relates to recent Forster reso¬ 
nance energy transfer (FRET) studies on folded and intrinsically disordered proteins 
(IDPs) and can account for the effects of solvent quality as well as the denaturant 
concentration on the reconfiguration dynamics. Following an ansatz that relates the 
strength of the harmonic confinement (kc) with the internal friction of the chain 
SDCRIF can convincingly reproduce the experimental data and explain how 
the denaturant can change the time scale for the internal friction. It can also predict 
near zero internal friction in case of IDPs. In addition, our calculations show that 
the looping time as well as the reconfiguration time scales with the chain length N 
as ~ where a depends weakly on the internal friction but has rather stronger 
dependence on the solvent quality. In absence of any internal friction, a = 21/ 1 

and it goes down in presence of internal friction, but looping slows down in general. 
On the contrary, poorer the solvent, faster the chain reconfigures and forms loop, 
even though one expects high internal friction in the collapsed state. However, if the 
internal friction is too high then the looping and reconfiguration dynamics become 
slow even in poor solvent. 
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I. INTRODUCTION 

Among the polymer rheologists the notion of internal friction associated with a single 
polymer chain is more than twenty hve years old UHl. Surprisingly it is only very recently 
that the topic has gained attention in the chemical and biophysics community. Although 
some earlier works in chemical physics community did point out the importance of internal 
friction in single chain dynamics nu but did not receive the attention it should have. 
However recent experiments on looping and folding dynamics [7HI2] in proteins have indi¬ 
cated the non-negligible role of internal friction. In this context it is worth mentioning that 
loop formation between any two parts of a bio-polymer [T3H21] is supposedly the primary 
step of protein folding, DNA cyclization and since internal friction affects looping it is obvi¬ 
ous that measurements of folding rates in proteins would predict the importance of solvent 
independent internal friction as well [25]. Subsequently not only these recent experiments 
[TTl [121126| on polypeptides and proteins showed internal friction to play a pivotal role in 
the dynamics but also motivated theoretical chemical physicists to come up with statistical 
mechanical models for single polymer chain with the inclusion of internal friction [2ZH3I] 
and apply these models to investigate the loop formation dynamics. Other than model build 
up there have been attempts to elucidate the origin of internal friction in proteins based on 
computer simulation studies [32H38] . A careful literature survey would reveal that it was de 
Gennes [32], who introduced the concept of internal viscosity at the single chain level. Rabin 
and Ottinger proposed an expression for the relaxation time, r^ei associated with internal 
viscosity following an idea of de Gennes [39], which is Trei = /ksT^rfs+rji) where, R = aN'^ 
and a, N are the monomer size and chain the length respectively, v is the Flory exponent 
[30Hl2]. Therefore in the limit solvent viscosity 0, it has a non-zero intercept propor¬ 
tional to the internal viscosity rji. This is what exactly seen in recent experiments where the 
plot of reconhguration time vs solvent viscosity has a hnite intercept equal to the time scale 
for the internal friction. To the best of our knowledge so far all the theoretical attempts 
on loop formation in a single chain with internal friction have been restricted to 9 solvent, 
with v = 1/2, when the chain behaves ideally. But the experimental conditions remain 
close to a good solvent rather than a 9 solvent. There have been few theoretical studies to 
elucidate the effect of solvent quality on loop formation in single polymer chain jlSl - HS] and 
unfortunately there exists almost no theoretical study to analyze the combined effect of sol- 
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vent quality and internal friction on the ring closure dynamics in polymer chains other than 
the very recent simulation by Yu and Luo |1U]. Apart from the solvent quality, denaturant 
concentration does play an important role in ring closure dynamics [211 SZl of proteins as it 
profoundly affects the compactness of the protein and routinely used in experiments. Very 
recently Samanta and Chakrabarti EB used a compacted Rouse chain model with internal 
friction to infer the role of denaturant on ring closure dynamics. But the model works only 
in the 6 solvent condition where a phantom Rouse chain description works. Experiments 
have been performed with proteins away from the 6 conditions where excluded volume in¬ 
teractions along with the internal friction play important roles. To take care of excluded 
volume interactions one has to go beyond phantom chain description but then the many 
body nature of the problem does not allow an analytical solution. One possibility would be 
to perform computer simulation as is recently done by Yu and Luo |16]. Other possibility 
would be to work with a polymer chain where the excluded volume interactions are taken 
care of at the mean held level. We take the second route and propose a very general model 
that takes care of solvent quality as well as denaturant in addition to internal friction. We 
call it “Solvent Dependent Compacted Rouse with Internal Friction (SDCRIF)”. To take 
care of solvent effect we closely follow the work of Panja and Barkema jlH] where an ap¬ 
proximate analytical expression for the end to end vector correlation function for a hexible 
chain in an arbitrary solvent was proposed based on a series of computer simulations. The 
expression carries a parameter u similar in the spirit of Flory exponent [201 - 12211201100] which 
takes care of solvent quality. A value of z/ = 1/2 corresponds to a 0 solvent and in that case 
the correlation function is exact and reduces to the text book expression for the ideal chain 
[SHEg, on the other hand u = 3/5 (0.588 more precisely) [28] corresponds to a self avoid¬ 
ing hexible chain (good solvent) as is the case with real polymers. Importantly the same 
expression can be used for a range of values of z/ corresponding to diherent solvent qualities. 
Very recently such a mean held Flory exponent based model has been used to describe ring 
polymer dynamics [53] . Next is the inclusion of internal friction which is done similarly as in 
case of a phantom polymer chain [52] . Further to take care of the denaturant which controls 
the compactness of the protein a conhning harmonic potential with force constant kc is used 
as is done in a very recent study by the authors [SI] and also in the context of dihusing 
polymers in microconhnements [22] or in case of bubble formation in double stranded DNA 
[56] . So the novelty of SDCRIF remains in its applicability for a range of solvent quality, 
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FIG. 1: Solvent Dependent Compacted Rouse with Internal Friction (SDCRIF) 

denaturant concentration and in addition it also takes care of internal friction {^int) when 
reqnired. SDCRIF rednces to Ronse model in the limit u = 1/2, kc = = 0 and to “ 

Compacted Rouse with Internal Friction (CRIF)” [30] in the limit kc ^ 7^ 0, z/ = 1/2. 

The looping dynamics is studied within Wilemski Fixman (WF) framework [T31107| with 
SDCRIF, assuming the polymer chain to be Gaussian. It is worth mentioning that WF 
formalism has extensively been used to calculate the average ring closure or looping time 
in presence of hydrodynamic interactions by Chakrabarti [58] and to elucidate the effect of 
viscoelastic solvent [59l [60] and even applied to ring closing of a semiflexible chain [61] . 

The arrangement of the paper is as follows. The model SDCRIF is introduced in section 
II. Section III deals with the method used to calculate the reconhguration and ring closure 
time. Results and discussions are presented in section IV and section V concludes the 


paper. 










































II. SOLVENT QUALITY DEPENDENT COMPACTED ROUSE WITH 
INTERNAL FRICTION (SDCRIF) 
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In the Rouse model, the polymer chain is treated as a phantom chain [SB [52], where 
the hydrodynamics interactions and the excluded volume effects are not present. If Rn{t) 
denotes the position of the monomer at time t of such a chain, made of (iV+1) monomers, 
n varying from 0 to N. The equation describing the dynamics of the chain is the following 


cdRnit) _ d‘^Rn{t) 


( 1 ) 


Where ^ denotes the friction coefficient and k = is the spring constant where b is the 
Kuhn length. f{n,t) is the random force with moments 


(/(n,t)) = 0 ,{fo,{n,ti)f^{m,t 2 )) = 2^kBT6a^6{n - - ta) (2) 

The above equation is solved by decomposing it into normal modes as follows R„(f) = 

OO 

Xo + 2 -^p(^)cos(^). In the normal mode description the above force balance equation 

p=i 

transforms to the following 


= -KX,{t) + /,(() (3) 

The time correlation function of the normal modes is 

{Xpo,{Q)Xqp{t)) = ^^5p/a/3exp(-f/rp) (4) 

fvp 

The above expression is for a phantom polymer chain, a chain in 6 solvent. In reality 6 
condition is rarely achieved and the chain behaves as in a good solvent where the chain 
swells due to the excluded volume interactions. For example Buscaglia et. al showed 071 
that to a 11 residue polypeptide chain, 6M GdmCl and 8 M urea are good solvents but 
aqueous buffer is very close to a 6* solvent. It is also not uncommon to encounter a situation 
representing more of a bad solvent where the chain has a collapsed conformation. To take 
care of solvent quality we introduce an exponent u similar in the spirit of Flory exponent 
isa in the above expression for the time correlation function of the normal modes. This 
inclusion closely follows the work of Panja and Barkema jUj. Apart from the solvent quality 
the presence of denaturant also alters protein conformation. Recently such effects have been 
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Parameter 

SDR 

SDRIF 

SDCR 

SDCRIF 


^SDR ^ 2Ni 

tSDRIF _ cSDR 1 27r2p2-'+l5,„t 

cSDCR _ cSDR. 

Sp Sp 

cSDCRIF _ cSDRIF 

Sp Sp 

kp 

uSDR _ 6w^kBTp^‘'+^ 

Fp — 

j^SDBJF _ j^SDR 

^SDCR = 2Nkc + kSDR 

^SDCRIF ^ 2Nk, + kl^R 

Tp 


^SDRIF 1 SDRIF 

.SDCR 

.SDCRIF , .SDCRIF 

P pM- + l 

~p ' ~int 

'p 

'p J^SDCRIF 1 ‘int 

'^int 

0 

^SDRIF _ ^int 
‘int k 

0 

.SDCRIF _ iint 

Cnt k+kcN'-“'+^ 


TABLE I: List of parameters for SDR, SDRIF, SDCR and SDCRIF 


taken care of by introducing a confining potential — 0)^) where, kc is the spring 

constant. In addition such a chain can posses what is known as the internal friction, ^int and 
can be termed as “Solvent Dependent Compacted Rouse with Internal Friction (SDCRIF)”. 
Earlier works discuss in detail how the inclusion of internal friction is done [29l [5l] • Like the 
Rouse model the end to end distribution of the polymer chain remains Gaussian in SDCRIF 
as well. Experiments also showed that Gaussian Distributions of end to end distances of 
proteins are not too bad approximations [IIIE2]. Any deviations from the Gaussian behavior 
is negligibly small |63] . For SDGRIF the basic structure of the correlation function remains 
the same but has three extra parameters z/, ^int and kc- 

{Xpa{0)Xq^{t)) = (5) 

Where, + 2Nk, and = 2N^ + The relaxation 

cS DC RIF 

time for mode is = ^sdcrif ■ As expected with v = 1/2 which corresponds to 6 

solvent, SDGRIF gives back “compacted Rouse with internal friction (GRIF)” |3T], when kc 
and ^int both are nonzero. A more realistic situation would be z/ = 3/5 corresponding to good 
solvent. If the confining potential is removed, the model no longer represents a compacted 
chain but still remains solvent dependent Rouse and can be termed as “solvent dependent 
Rouse with internal friction (SDRIF)”. Further if the internal friction S,int is ignored but 
v ^ 1/2, the model reduces to solvent dependent Rouse (SDR), and as expected, in the 
limit when kc both are zero and z/ = 1/2, Rouse chain is recovered. Another model 
would be “solvent dependent compacted Rouse (SDGR)” when ^int = 0, but kc ^ 0. SDGR 
can account for denaturant effect even in the absence of internal friction which we have 
used later to reproduce the reconfiguration time of an IDP in high denaturant concentration 
when internal friction should be negligible. Fig. ([^ summarises the models and Table. 
depicts all the models along with the parameters describing the models. 
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III. CALCULATION METHODS 


A. Reconfiguration time 


The time correlation function of the end-to-end vector is calculated from the correlation 
of normal modes as follows 


Therefore, 


/ TD D /i\\SDCRIF -| n Shj^T . , 

p=i > 


SDCRIF\ 
P ' 


( 6 ) 


0Aro(O) — (^Rno^ 


SDCRIF 


eq 


= 16E 


SksT 


USDCRIF 

=1 N 


(7) 


The exact expression for is not analytically trackable. However for 9 solvent 

when z/ = 1 / 2 , it has an analytical expression [3l] 


2 \CRIF 26a/3/csT ^ ^ Nby/¥c 

^ 7 ^— tClTZh^~ 


•^^ 0 /eg,.= 1/2 


Vk~c 


'2^3^' 


( 8 ) 


Fortunately for SDRIF, (-R^o)eg analytical expression too 


o \ SDRIF 


■)2—2u 


(-1 + 2i+2^)C[l + 2 z/]62a2 


TT^ 


(9) 


Where, C[x] = / ^ 7 —j- is Riemann Zeta function and F(a;) is the gamma function. 

Reconfiguration time tno is obtained by integrating the normalized (j)m(t) 


Tm — J dt(j)NQ{t) 
0 


( 10 ) 


Where, 


B. Looping time 


Wilemski Fixman (WF) approach is a widely accepted method to calculate the time 
required to form a loop between two parts of a Gaussian polymer chain [13]. This formalism 
gives the following expression for the looping time between two ends of the chain. 


loop _ 

Wo ~ 



/ GAro(t) 
\GAro(c)o) 



( 11 ) 












Where, Cjsioit) is the sink-sink correlation function given by 


Cmit) — / dRm / dRNo,oS{RNo)G{RNo,t\Rm,o,^)S{Rm,o)P{dij^Qfi) 


( 12 ) 


G{Rm, tlRmfi, 0) is the conditional probability of the polymer to have end-to-end distance 
Rno at time t, which was Rno,o at time t = 0. 


\ 3/2 


G{Rm,t\RNo,o,0) — 


o /d2 \SDCRIF 

.^^{Pm)eq j 




exp 


^{Rm — 4>No(t)Rm,o) 

, SDCRIF 


2Wo)e7 

(13) 

When the sink function S{Rmn) [61H66] is chosen to be a delta function then the looping 
time of SDCRIF has the following expression |S7| 


T 


loop,SDCRIF _ 
NO ~ 


dt 


exp[-2xo^mit)/i^ - ^%Q{t))]sinh[{2xo^No{t ))/(1 - 077 )] 


(2xo07vo(t))Vl - 07(^) 


Where, 


Xo = 


3a^ 


2 (-R’no) 


eq 


(15) 


IV. RESULTS 


A. Equilibrium end to end distribution 


The equilibrium distribution of the vector connecting end-to-end monomers of a Gaussian 

3/2 


chain is given by P{R]sfo) = 


2n 


{^%o) 


exp 


3^^o 

^Ko) 


, where, {R%q) denotes the average 


equilibrium end to end distance of the polymer. This expression holds for SDCRIF as 
well since the model is considered to be Gaussian in our description. In Fig. (|^ the 
equilibrium distribution of the end-to-end distance is shown for the polymer at different 
solvent quality in absence of the conhning potential. As expected the distribution plot is 
broader for the polymer in good solvent (z/ = 3/5) in comparison to the polymer in poor 
solvent (z/ = 1/3). From good to 6 to bad solvent the most probable end to end distance 
or the peak position shits to lower value as a signature of swelled to collapsed transition. 
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The parameters have been chosen in accordance to the work by Schuler’s group m on 
cold shock protein (Csp) and prothymosin a (ProTa) and has been mentioned in all the 
hgures. In Fig. (|^ the same distribution is shown for good solvent but in presence of three 
different values of kc. As the kc controls the compactness of the polymer, increasing the 
value of kc should result in the higher degree of internal friction. Since there is no hrst 
principle relation between kc and we have used an ansatz in our calculation [IHEI]. 
The ansatz connects kc with ^int as follows kc = kc + kcfl{cQ + Cirih + C 2 n^ + ....), where, 
kc,o = and = 6nt,o(co + CiUfe + C 2 nl + ....). is the zeroth order approximation 

to internal friction which only accounts for the interactions between a monomer with its two 
nearest neighbours, Tint,o is the corresponding time scale, rib is the number of non-adjacent 
monomers contributing to the internal friction and A is considered typically in the order of 
Ti ‘^whereas other parameters are constants with no or very weak dependence on chain 
properties. Therefore kc has a part proportional to However, in the limit ^intfl = 0, 

kc = kc- If the value of kc or the strength of the conhnement increases the end to end distance 
probability distribution becomes narrower, a sign of more compacted polymer chain. This 
can clearly be seen from the plot of end to end distance distribution at different values of kc 
depicted in Fig. ([^. Thus kc and u plays signihcant roles in deciding the width and height of 
the end to end distribution. But the the parameters have different physical signihcance and 
play different roles in the chain dynamics. The exponent z/ takes care of the quality of the 
solvent around the chain and has no dependence on ^int and no or very weak dependence 
on denaturant concentration m- But kc on the other hand takes care the effect of the 
denaturant on the chain and has ^int dependence as well. Here we also attempt to make a 
qualitative comparison of our model SDCRIF with the very recent simulation from Luo’s 
[l6] group. Fig. 3. of [16], shows similar trend of the end to end probability distribution on 
changing the parameter A, a measure of attractive interaction between the beads. Higher 
the value of A poorer the solvent, narrower the distribution. 

B. Relaxation time for the Normal modes 

In order to calculate reconhguration time one has to calculate the correlation function 
between normal modes dehned in Eq. The functional from of is given in Table. 

The lower normal modes contribute largely to looping and reconhguration dynamics [57]. A 
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FIG. 2: End-to-end distribution at different solvent quality with fee = 0. The values of parameters 
used are N = 66, b = 3.8 x ^ = 9.42 x ks = 1.38 x JK~^ and T = 300iir. 

plot of vs z/ for the first normal mode is shown in Fig. Q. This shows how 

changes with the solvent quality if all the other parameters remain fixed. The plot shows 
that (for p = 1) increases with u. Thus for poor solvents corresponding to lower 

value of u, is lower and as one approaches good solvent increases resulting 

in slower relaxation. This would definitely lead to slower reconhguration and looping time 
in good solvent as compared to bad or poor solvent. Similar trend for the third normal 
mode can be seen from the plot in the inset of the Fig. (|^. Recent simulation also confirm 
this |in] and shows an order of magnitude increase in the looping time in the good solvent 
as compared to the poor solvent. The reconfiguration times of a chain for different values 
of V are shown in Fig. (|^. It can be seen that at a value p/po = 1, fhe reconfiguration time 
increases by a factor of 6 on changing z/ = l/3toz/ = 3/5. 

C. Solvent viscosity dependence 

One aspect of the reconfiguration time is its dependence on solvent viscosity [29]. Experi¬ 
ments show a linear dependence of reconhguration time on solvent viscosity HI]. In Fig. (g 
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FIG. 3: End-to-end distribution at different values of kc for good solvent {u = 3/5). The values 
of parameters used are N = Q6, kc = 0, b = 3.8 x ^ = 9.42 x 10~^'^kgs~^ , ks = 

1.38 X 10"23 jK-i and T = 300K. 

the reconfiguration times are plotted against the solvent viscosity at different values 

of u in absence of any conhning potential {kc = 0) and non zero internal friction {^int ^0). 
The hgure clearly shows increases linearly with increasing normalized solvent vis¬ 
cosity r]/riQ, with a positive intercept in the limit —)■ 0, rjo being the viscosity of water. 

What is surprising is the value of intercept is practically independent of solvent quality. This 
intercept corresponds to the internal friction which is completely dry in this case. The trend 
can be htted with an almost analytically exact expression + C{h')T^^^, where 

^SDR _ and C'(z/)s have v dependence and are 0.73, 0.82 and 0.87 for the poor, 

9 and good solvent respectively. Thus poorer the solvent lower the slope. In poor solvent 
reconhguration time is not only fast but also less affected by the viscosity of the solvent. 
This is presumably because of the collapsed form of the polymer chain. Fig. (|^ is the plot of 
vs 77/770 in presence of a hxed conhning potential /Cc,o but at three different values of 
z/ covering a range of solvent qualities. Interestingly in this case the intercepts are found to 
be different when the intercepts were magnihed as can be seen in the inset of Fig. ([^. This 
is because has kc as well as v dependence as can be seen from the Table. |lj So even 
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V 


FIG. 4: vs V for p = 1 and p = 3 (inset). The values of parameters used are N = 66, 

kc = kc,o, kc = 0,b = 3.8 x ^ = 9.42 x 10-^^kgs-\ ^int,o = 100 x^,kB = 1.38 x IQ-^^ JR-^ 

and T = 300iir. 

with the same value of k^, changing u would result a change in the effective force constant, 
keff — k + kcN ‘^'^^^of the chain and hence a change in the value of the intercept. In 
all the cases ^int remains the same and the time scales associated with the internal friction 
are changed by very small amount which means the dependence of on u is 

very weak and it cannot reproduce the denaturant effect observed in the experiments, which 
is why we need k^ and the ansatz which connects kc with ^int to take care of the denaturant 
effect on a polymer without changing the solvent quality. This can be seen in Fig. ([^ where 
we have looked into the viscosity dependence of in different values of kc for good 

solvent {u = 3/5). In this case the change in the values of the intercepts is very evident 
because of the different degrees of compactness of the chain. Another important observation 
is the case with higher kc value has a higher intercept but it is less steeper making the 
change in the reconhguration time less prone to the viscosity of the solvent. This is how 
kc can replicate the experimental results m where it has been found that if denaturant is 
introduced to a polymer solution, the polymer becomes less compact which results in lower 
value of internal friction ^int- kc controls the compactness of the protein thus connects to the 
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FIG. 5: Reconfiguration time vs solvent viscosity at different solvent quality. The values 

of parameters used are N = 66, kc = 0, b = 3.8 x ^ = 9.42 x 10~^'^kgs~^, ^intfl = 100 x 

kB = 1.38 X JK-^ and T = 300R:. 

denaturant concentration. The higher the denaturant concentration lower the kc value. This 
is also reflected in our ansatz as increasing denaturant concentration would mean smaller 
Ub and lower internal friction But kc does not speak for the solvent quality rather it is 
u which accounts for that. In the later section we used SDCRIF to compare experimental 
data on cold shock protein and prothymosin a (ProTa) [TT] and have confirmed this. We 
have also looked at the viscosity dependence of the looping time and have found it to be 
~ r]^ with (3 <2 [29l EH] . A detailed study on the viscosity dependence of the looping time 
can be found in one of our earlier works |29j . 


D. Chain length dependence 

Fig. (|^ is the log-log plot of the reconhguration time for a chain with kc = = 0 

vs chain length N of the polymer for different values of z/, the parameter accounting for the 
solvent quality. There is a general trend, ~ as This result is 

expected and can be predicted by looking at the scaling of the end to end vector correlation 
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FIG. 6: Reconfiguration time vs solvent viscosity at different solvent quality. The 

values of parameters used are N = 66, kc = kcfl, kc = 0, b = 3.8 x ^ = 9.42 x 10~^‘^kgs~^, 

O = 100 x^,kB = 1.38 X and T = 300R:. 


function with the chain length. When the same calculations were done in presence of internal 
friction as shown in Fig. (|^ the dependence on the chain length becomes weaker in general 
but the relative trend remains the same, in poor solvent dynamics is faster. Although the 
chain reconhguration in general become slower due to the presence of internal friction. To see 
how the chain dynamics changes in presence of the conhning potential the same calculations 
were again performed for three different values of kc but putting internal friction ^int = 0 


in good solvent = 3/5) as shown in Fig. (10). Here as ^int = 0 therefore, /cc,o = 0 
and kc = kc whereas kc is scaled as /cc,o- This result is very interesting as it can be seen 


from the Fig. (10), that the N dependence remains practically unchanged. The reason for 
this as follows. We restrict ourselves to a value of kc such that the effective force constant 
kcff — k + kcN ‘^^~^^/is in the same order of magnitude of k. For example a choice 
of kc — gives fee// — k{l + Now since lower normal modes contribute 

mostly in reconhguration dynamics 153 , taking the contribution from p = 1 normal mode 
gives keff — 2k and thus leads to unchanged scaling of the reconhguration time with the 
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FIG. 7: Reconfiguration time vs solvent viscosity at different values of kc for good 

solvent (z^ = 3/5). The values of parameters used are N = 66, kc = 0, b = 3.8 x 
^ = 9.42 X IQ-'^'^kgs-^, ks = 1.38 x and T = 30077. 


chain length N for a given value of v. Similar calculations have been done for looping time 
based on the WF approximation na as can be seen from Fig. © and Fig. ( |l^ . The 
same arguments hold for the looping time and the trends are similar. Scaling of the looping 
time as N'^ for the Rouse chain [58l |69] as obtained earlier can be conhrmed with a choice 
of z/ = 1/2, fee = 0 and ^int = 0 and is shown in Fig. ©■ 


E. Comparison with experiments 

In this section we use SDCRIF to compute the reconhguration time and compare it 
with the one measured in the recent Forster resonance energy transfer (FRET), nanosecond 
fluorescence correlation spectroscopy and microfluidic mixing based study on cold shock 
protein from Thermotoga maritima (Csp) labeled at positions 2 and 68 with Alexa 488 and 
Alexa 594 as donor and acceptor, respectively [IIIITD]. Thus it can be approximated that the 
labelling are at the two ends of the protein. Their study revealed important role of internal 
friction in unfolded small cold shock protein and conhrmed a time scale of about ~ 5 — 50 
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FIG. 8: Log{T^Q^) vs Log{N) at different solvent quality. The values of parameters used are 
= 0, 6 = 3.8 x ^ = 9.42 x lQ-^‘^kgs-^, ks = 1.38 x 10-‘^^JK-\ = 0 and T = 300^. 


ns for the internal friction. Moreover they found the internal friction to strongly depend 
on the denaturant concentration. Nuclear magnetic resonance and laser photolysis methods 
also have confirmed the effects of denaturants by showing the rate of intrachain contact 
formation in unfolded state of carbonmonoxide-liganded cytochrome c (cyt-CO) to increase 
with increase in denaturant concentration na. Higher the denaturant concentration lower 
the internal friction. On the other hand for an intrinsically disordered protein (IDP) such 
as C-terminal segment of human prothymosin a (ProTa) magnitude of internal friction 
is smaller [m El and at high denaturant concentration it is negligibly small. This is 
presumably due to exposed hydrophilic and charged residues resulting expansion of the 
protein. But in native buffer the time scale for internal friction is ~ 6 ns which on addition 
of excess salt like KCl shows ~ 3 times increase in internal friction due to the collapse of 
the IDP. For both the cases the reconhguration time measured has a linear dependence on 


solvent viscosity with an intercept equal to the time scale for the internal friction. Fig. (13) 


and Fig. (14) show the experimental data and our calculation based on SDCRIF of Csp and 


ProTa respectively. All the calculations are performed with u = 3/5, value corresponding 
to good solvent, which is particularly valid for IDPs. On the other hand, the parameter 
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FIG. 9: Log{T^Q^^^) vs Log{N) at different solvent quality. The values of parameters used are 
kc = 0, b = 3.8 X ^ = 9.42 x iintg = 100 x ks = 1.38 x 10“^^Ji7“^ and 

T = 30077. 

kc in our model takes care of changes associated with the denaturant concentration and 
^int is the corresponding internal friction which is considered to be zero for the ProTa in 
high denaturant condition (OMGdmCl) when time scale associated with internal friction was 
found to be negligible experimentally. In this case kc = kc and kc is scaled as 5.5/cc,o- Table. 
0and Table. ^ depict the parameters used in our calculation and show the theoretically 
calculated internal friction time scales are in excellent quantitative agreement with that 
measured experimentally. 


V. CONCLUSIONS 

Motivated by recent experiments m on cold shock proteins and intrinsically disordered 
proteins (IDPs) we have analyzed the effects of denaturant and the solvent quality on the 
reconhguration and looping dynamics of a chain with internal friction by using an extended 
Rouse chain model with internal friction. The model termed as “Solvent Dependent Com¬ 
pacted Rouse Chain (SDCRIF)” takes care of solvent quality through a Flory type exponent 
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FIG. 10: vs Log{N) at different values of kc for good solvent {u = 3/5). The values 

of parameters used are b = 3.8 x ^ = 9.42 x 10~^^kgs~^, iintfl = 0, fes = 1.38 x 10“^^ 

and T = 300111. 

V and the effects of denaturant concentration are taken care by the strength of a harmonic 
conhnement kc of the chain. Following an ansatz we farther relate kc with the internal 
friction. This assnres a non zero intercept in the plot of reconhgnration time vs solvent 
viscosity as fonnd in experiments dUEo] and also makes it denaturant concentration de¬ 
pendent. Here we would like to point out that mere “Compacted Rouse with Internal Friction 
(CRIF)” [m El] can not convincingly account for the changes in reconhgnration time due 
to change in solvent quality and for this we need a parameter z/. Also the parameter kc 
should be coupled strongly with the internal friction ^int and this is done through the ansatz 
mentioned in the result section. For folded protein like cold shock protein the magnitude of 
the internal friction is high and our theory also reproduce this with a choice of high value 
of kc and ^int- Values of these parameters also change on changing the denaturant concen¬ 
tration. Typically the value of the number of monomers contributing to internal friction, 
rif) also increases as the denaturant concentration decreases and the protein collapses which 
demands a higher magnitude of kc to be used to reproduce the experimental data. On the 
other hand for the intrinsically disordered protein prothymosin a (ProTa) magnitude of 
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FIG. 11: Log{Tj^Q^’^^^) vs Log{N) at different solvent quality. The values of parameters used 
are kc = 0, b = 3.8 x 10~^^m, ^ = 9.42 x 10“^^A:g's“^, ^int,o = 0, fes = 1.38 x and 

T = 300ii:. 


kc and the internal friction ^int are low. Only in presence of a salt the protein collapses 
and then a higher magnitude of kc can reproduce the experimental data. This highlights 
the novelty of our model which is applicable to a wide range of denaturant concentration, 
solvent quality and protein types. Whereas kc in our model effectively renormalizes the force 
constant of the chain, u on the other hand takes care of the quality of the solvent around. 
Both the parameters are essential to explain the experimental data on reconhguration times 
of proteins, so as the azsatz that relates kc with the internal friction ^int of the chain. An¬ 
other issue is the relative dependence of the reconhguration and looping time on and v. 
For the range of values of ^int we have used throughout our calculation the reconhguration 
time (so as the looping time, not shown) depends only weakly on internal friction and 


rather strongly on the solvent quality, i.e. v. This can be seen in Fig. (15) (where c{v) 
is taken to be 0.80 for simplicity) but on increasing the internal friction magnitude by two 
orders of magnitude the change due internal friction become visible. For example if the ^int 
is taken to be 3000^ for poor solvent (z/ = 1/3) the time scale for internal friction become 
(~ 340 ns) which is ~ 3 times higher than the same (~ 100 ns) in good solvent (z/ = 3/5) for 
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FIG. 12: vs Log{N) at different solvent quality. The values of parameters used are 

kc = 0, b = 3.8 X ^ = 9.42 x ^mt,o = 100 x ks = 1.38 x 10“^^and 

T = 300ii:. 

6nt = 100^. This is in the same spirit as seen in recent simulation |1B], where the looping 
time passes through a minima while plotted against the parameter A, a combined measure 
of solvent quality and internal friction. 

We would like to conclude by pointing out that in reality internal friction has contri¬ 
butions from hydrogen bonding, other weak forces and specially torsion angle rotations in 
proteins and have already been investigated in atomistic simulations [35flS7] . Taking these 
contributions explicitly beyond the scope of this study. However it would be worth incor¬ 
porating internal friction in a model of polymer with torsion and semiflexibility [52] and 
explore the physics involved. Work along this direction is under progress. 


VI. ACKNOWLEDGEMENT 

Authors thank IRCC IIT Bombay (Project Code: 12IRCCSG046), DST (Project No. 
SB/SI/PC-55/2013) and CSIR (Project No. 01(2781)/14/EMR-II) for funding generously. 







21 




(a) (b) 




FIG. 13: Comparison between theoretical and experimental reconfiguration time vs solvent vis¬ 
cosity data for the cold shock protein (Csp). (a) 1.3M GdmCl, (b) 2M GdmCl, (c) 4M GdmCl 
and (d) 6M GdmCl. The values of parameters used are N = 66, kc = 0, b = 3.8 x 
^ = 9.42 X IQ-^'^kgs-^, Cmi,o = 100 x ks = 1.38 x and T = 300K. 


Denaturant concentration 

he 

^int 

Tint (Theory) 

Tint (Experiment) 

1.3M GdmCl 

3.0kc,o 

9*0<^2nt,0 

~ 47 ns 

~ 42 ns 

2.0M GdmCl 

2.5/cc,o 


~ 27 ns 

~ 25 ns 

4.0M GdmCl 

2.0/cc,o 

3.0i^j,^t^o 

~ 17 ns 

~ 12 ns 

6.0M GdmCl 

l.bkcfi 


~ 9 ns 

~ 5 ns 


TABLE II: Comparison with experimental data on cold shock protein: The values of parameters 
used are N = 66, kc = 0, b = 3.8 x 10“^°m, ^ = 9.42 x 10~^^kgs~^, iintfl = 100 x ks = 
1.38 X 10-23 JiL-i and T = 300A:. 
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FIG. 14: Comparison between theoretical and experimental reconfiguration time vs solvent vis¬ 
cosity data for IDP ProTa. (a) native buffer, (b) IM KCl and (b) 6M GdmCl. The values 
of parameters used are = 110, 6 = 3.8 x ^ = 9.42 x , iintg = 100 x 

ks = 1.38 X 10"23jic-i and T = 300iP. 
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kc 

^int 
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Native buffer 

kc — 0? 6/cc,o 
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IM KCl 

kc — 0? 8/cc,o 


~ 12 ns 
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0 

0 
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FIG. 15: 3D Plot of reconfiguration time vs v and ^int- (a) ^int upto 10^ x (b) ^int upto 10^ x 
The values of parameters used are N = 66, kc = 0, b = 3.8 x ^ = 9.42 x 10~^‘^kgs~^, 

ks = 1.38 X and T = 300iF. 
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